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We construct an energy-momentum tensor on the lattice which satisfies the appropriate Ward 
Identities (WIs) and has the right trace anomaly in the continuum limit. It is defined by imposing 
suitable WIs associated to the Poincare invariance of the continuum theory. These relations come 
forth when the length of the box in the temporal direction is finite, and they take a particularly simple 
form if the coordinate and the periodicity axes are not aligned. We implement the method for the 
SU(3) Yang-Mills theory discretized with the standard Wilson action in presence of shifted boundary 
conditions in the (short) temporal direction. By carrying out extensive numerical simulations, the 
renormalization constants of the traceless components of the tensor are determined with a precision 
of roughly half a percent for values of the bare coupling constant in the range 0 < < 1. 


I. INTRODUCTION 

On the lattice the Poincare group is explicitly broken 
into discrete subgroups, and the full symmetry is recov¬ 
ered only in the continuum limit. As a consequence, a 
given definition of the energy-momentum tensor needs 
to be properly renormalized to guarantee that the asso¬ 
ciated charges are the generators of the Poincare group 
in the continuum limit, and that the trace anomaly is 
correctly reproduced. 

In order to construct the renormalized energy- 
momentum tensor, the way to proceed is to impose suit¬ 
able WIs at fixed lattice spacing that hold up to cutoff 
effects which vanish in the continuum limit [l[. This pro¬ 
gram can be realized in practice if the WIs involve corre¬ 
lators which in turn are simple enough to be computed 
numerically with good precision. 

When the theory is considered in a finite box, the Eu¬ 
clidean Lorentz symmetry is also softly broken by its 
shape. If the length in one (temporal) direction Lq is 
chosen to be shorter than the typical scale of the theory 
(thermal theory), interesting WIs follows They 

become particularly simple and of practical use if the pe¬ 
riodicity axes are tilted with respect to the lattice grid, 
i.e. if the hard breaking of the Poincare symmetry due 
to the lattice discretization and the soft one due to the 
finite temporal direction are not aligned. This set-up has 
a natural implementation in the Euclidean path-integral 
formulation in terms of shifted boundary conditions @0 . 

Here we define the renormalized energy-momentum 
tensor of the Yang-Mills theory non-perturbatively by 
working in this framework. This is achieved by supple¬ 
menting the set of WIs found in Refs a a with a new 
one which guarantees that the correct trace anomaly is 
reproduced in the continuum limit. 

We implement this strategy for the SU(3) Yang-Mills 
theory regularized with the standard Wilson action. We 


carry out extensive numerical simulations, and compute 
with high precision the renormalization constants of the 
traceless components of the tensor. The numerical de¬ 
termination of the renormalization constant of the trace 
part requires additional simulations, and it is left for a 
forthcoming publication. Over the last year an alterna¬ 
tive method, based on the Yang-Mills gradient flow, has 
also been proposed for renormalizing non-perturbatively 
the energy-momentum tensor J6 h8J. 


II. WARD IDENTITIES IN PRESENCE OF A 
NON-ZERO SHIFT 

In this section we consider the SU(3) Yang-Mills the¬ 
ory in the continuum. The definitions of the action and 
of the partition function Z are reported in Appendices [A] 
and |B] together with other conventions. Here we are in¬ 
terested in the thermal theory defined in the path integral 
formalism with shifted boundary conditions 

A^(L 0l x) = A M (0,a: - L 0 £) (1) 

along the compact (temporal) direction of length Lq with 
shift ^ £ R 3 . The free-energy density is given by 

f ( L o, Z) = -j^ In Z(L 0 ,£), (2) 

where V is the spatial volume of the box. In the thermo¬ 
dynamic limit, which is always assumed in this section, 
the invariance of the dynamics under the SO(4) group 
implies @| 

f(L 0 ,t) = f(L 0 y/TTe,0). 


( 3 ) 
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When £ ^ 0, odd derivatives in do not vanish, and the 
following interesting relations hold M (® 0 7 ^ 0 ): 

Lo( T 0 k)z = —-gjr- lnZ(Lo,£) , 

— (O)£=L 0 (T 0 k(x 0 )O(0))t, c , (4) 

where T gl/ is the energy-momentum tensor, T = 
f g? 3 xT m „(x), O is a generic gauge invariant operator, and 
the subscript c indicates a connected correlation function. 
By deriving once with respect to Lq and to one obtains 
the relation (no summation over repeated k here) 0 

(T ofc ) € = {(T 00 ) c - (T kk )s } . (5) 

In the equations above and in the rest of this paper we 
focus on correlation functions of the energy-momentum 
tensor with gauge-invariant operators inserted at a 
physical distance. As reviewed in Appendix [UJ it is ap¬ 
propriate in those cases to consider the symmetric gauge- 
invariant definition of the energy-momentum tensor given 

by 


are two fields transforming as a two-index symmetric and 
a singlet irreducible representation of the SO(D) group 
respectively, and D = 4 — 2e. 

The field ris a dimension-four gauge invariant op¬ 
erator which is multiplicatively renormalizable. The WI 
in Eq. 0 fixes its renormalization constant to 1. This in 
turn implies that g q t^ v renormalizes as 

f 1 "I ^ 

\F° la F* a - -S^F^A = Z g {F- a F- a ( 12 ) 

1 A P a l 

a( 3 "aft | J 

where Z g is the renormalization constant of the coupling, 
see Appendix |D] By definingthe renormalization group 
invariant (RGI) operator as [9| 

2 (13> 

1 1 R 

~~pjfi/J.>'F a pF a p j ) 

we finally arrive to 


T? a l? a A T? a T? a 


RGI 


T„ 


-»{ 


T?d pa 

fia vol 






( 6 ) 


By deriving two times with respect to the shift compo¬ 
nents and by using Eq. 0 , one obtains 0 


(Tokh = ^'Z( T <* T °ikc 


5 _ & & 

£2 


(7) 


where on the r.h.s. the two fields are inserted at different 
times. Analogously one shows that (xo ^ 0) 

L 0 ( To fc (xo) T w ( 0 )) €>c = jb - ( T ofc)« + 

Ao^^(T ofc (xo)T ofc ( 0 )) € , c , ( 8 ) 
?fc 

which can also be put in the more suggestive form 


— ( T ) 1 9 
dtik { ^ (l + £ 2 ) 2 c >6 


2\3 


a+a 

6 


-{To k)t 


(9) 


A. Finiteness of T gv and trace anomaly 

In dimensional regularization, the energy-momentum 
tensor in Eq. 0 is decomposed as 


RGI 

KpKp) , ( 14 ) 

where bo is the first coefficient of the /3-function given in 
Eq. (ID41) . The field r is also dimension-four and gauge 
invariant, but it is a singlet under SO(D). Therefore it 
mixes with itself and with the identity operator. The 
Eq. 0 fixes the multiplicative renormalization constant 
to 1 , while a natural prescription for the identity sub¬ 
traction is 


— A jt/ — 2 &o 


pa pa 
flOL VOL 


~i s ' 


T R = T - (t )o (15) 


where (• • ■ )o indicates the vacuum expectation value for 
Lq —)> oo (zero temperature). This in turn implies that 
one can define 

{KpK 0 } R = Ze 1 {KpKp - {KpKp) o) , (is) 

and Eq. m becomes 


T 


R 


eZgZs 
2 D/j?*g 2 


{KpKp } 1 


(17) 


By using the result in Eq. (ID12I) of Appendix [0 one 
obtains the well known result for the trace anomaly [l(| 

El 


T gv — Tgv + 


where 


1 /IV 


So 


_ ) TT'a rpa r 771 a ttm 

O 1 11 ru ■T' i/rv t 1 0 IJ,V rv R -L , 


2 Dg { 


fia vol 


rpa T^a 
2 r aP r oLP 


D 


V-v-l afj- 1 - ol{l 


( 10 ) 


(ii) 


Tup = 4r R = -^3 {FZpFZp } R . (18) 

By defining the renormalization group invariant operator 
as 0 


f rpa jra 1 RGI _ 

aft-T ol/3 j 


bog ; 


r 771 a 771a 
t r OL(3 r a(3 


1 R 


(19) 


r = 
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we can finally write 

Tw = 4r R = ^{FZpFZp}™ 1 . (20) 

The WIs in Eqs. © and © fix unambiguously the renor¬ 
malization constants of the composite fields entering the 
energy-momentum tensor definition so that the correct 
trace anomaly is reproduced. The Eqs. (fTSl) (l20l) hold to 
all orders in perturbation theory. 


III. THE ENERGY MOMENTUM TENSOR ON 
THE LATTICE 


We regularize the SU(3) Yang-Mills theory on a finite 
four-dimensional lattice of spatial volume V = L 3 , tem¬ 
poral direction Lo, and spacing a. The gauge field sat¬ 
isfies periodic boundary conditions in the three spatial 
directions and shifted boundary conditions in the com¬ 
pact direction 


U li (L 0 ,x) = U ll (0,x-L o e), (21) 


where U^x o,%) are the link variables. The action is 
discretized through the standard Wilson plaquette 



EE 


1 - -!-Re Trl U, 


{u,Ax)} 


( 22 ) 


where the trace is over the color index, and /3 = 6/<?q 
with go being the bare coupling constant. The plaquette 
is defined as a function of the gauge links, and it given 

by 


U^{x) = U^{x) U v (x + afi) Ul(x + av) J7*( x) , (23) 

where fi, v = 0,...,3, /t is the unit vector along the di¬ 
rection g, and x is the space-time coordinate. The gluon 
field strength tensor is defined aSQ 


/L V 


(x) = ~ Tf{ [Q^{x) - Q vlt (x)] T a ) , (24) 


where 


Q,xv{x) = U^,(x) U u (x + afi) U^(x + av) Uf(x) 

+U v (x) Uj^(x — a/ i + av) Uf(x — afi) U^x — afi ) (25) 

+Uf(x — afi) Uf(x — afi — av) U )1 (x — afi — av) U„(x — av) 
+Uf{x — av) U /Ji (x — av) U u (x + afi — av) U'Ax) . 


1 We use the same notation for lattice and continuum quantities, 

since any ambiguity is resolved from the context. As usual, the 

continuum limit value of a renormalized lattice quantity, iden¬ 
tified with the subscript R, is the one to be identified with its 
continuum counterpart. 


The target energy-momentum tensor in the continuum 
is a gauge-invariant operator of dimension 4, which is 
a combination of a traceless two-index symmetric and a 
singlet irreducible representation of SO(4). When SO(4) 
is broken to the hypercubic group SWq, the traceless 
two-index symmetric representation splits into a sex¬ 
tet (non-diagonal components) and a triplet (diagonal 
traceless components). At finite lattice spacing, the 
energy-momentum tensor is thus a combination of gauge- 
invariant operators of dimension d, < 4 which, under the 
hypercubic group, transform as one of those two represen¬ 
tations and the singlet. In the SU(3) Yang-Mills theory 
there are only three such operators (no summation over 
repeated fi and v here) [3, El: 

T}^ = (i - ^,,)\{f^F-\ 
y o K J 

= ( 26 ) 

o 

r r\^\ £ ^ / tt'Cl r?a 1 77>a r?a \ 

fiv Vilv 2 1 ^[iol a Oifi-Ta(3 | 

y o k ^ J 

and the identity. The sextet and the triplet 

Tff renormalize multiplicatively, while the singlet Tff 
mixes also with the identity. The renormalized energy- 
momentum tensor can finally be written as 

T* = Z T {TW + z T Tf v + * s [T$ - (2$)o] } . (27) 

The renormalization constants Z T , z T and z s are finite 
and depend on only. At one loop in perturbation 
theory their expressions are EL El 

Z T {gl) = 1 + 0.27076 g 2 0 , 

z t («?o 2 ) = 1 -0.03008 gl, (28) 

( 2\ t>o 2 

z s(9o ) = y5o ■ 


A. Non-perturbative renormalization conditions 

The renormalization constants Z T , z T and z s can be 
determined non-perturbatively by requiring that on the 
lattice the WIs in Eqs. Q, ©, and © hold up to dis¬ 
cretization effects which vanish in the continuum limit. 
The renormalization constant of the sextet is fixed to 

be El ( see also El) 


Z t (9o) — ~ 


1 


A/ _ 

(Tok)t. 


(29) 


where the derivative in the shift in Eq. © is discretized 
by the symmetric finite difference 


A/ _ 1 ln r Z(£ 0 ,g-afc/£ 0 ) 

A4 2 aV n lz(L 0 ,€ + ak/L 0 ) 


( 30 ) 
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FIG. 1. Left: the derivative of Af/A^k (normalized to its Stefan-Boltzmann value) with respect to go as a function of the 
bare coupling. Right: (Tq^)^ normalized to its tree-level expression as a function of g$. The data are generated on lattices 
with Lo/a = 3 (blue), 4 (green) and 5 (red), L/a = 48, and £ = (1,0,0). Statistical errors are smaller than symbols. 


to ensure that discretization effects start at 0(a 2 ). In By choosing one possibility of discretizing Eq. the 

the thermodynamic limit, which is always assumed in singlet renormalization constant is fixed to be 

this section, the triplet is renormalized by requiring that 
Eq. (0 holds up to harmless discretization effects, i.e. 


(do ) 


i-4 2 (4 ] )$ 

a <!$>* - <2$> c ' 


(31) 


J 


(l+£ 2 ) 3 /rp[l]\ 

i' k ro fc /€ 

— 

(i-K 2 ) 3 / r W\ 

£' k VO k/£ 



£'=£+aic/ L 0 


£'=£—ak/L o 


~ s (l + £ 2 ) 2 (T m ) - -(T m \ - 

v s ’ \- L W-/£+ak/L 0 \- L W/£-ak/L 0 


At finite Lo, the renormalization constants depend on 
the bare coupling constant and on (a/Lo) 2 due to dis¬ 
cretization effects. Our prescription is to define them in 
the limi1@ Lq —> oo at fixed g$. 

IV. NUMERICAL COMPUTATION 

In this section we describe how the strategy outlined 
above has been implemented in practice to determine 
the renormalization constants Z T and z T . In all simu¬ 
lations the basic Monte Carlo step is a combination of 
heatbath and over-relaxation updates of the link vari¬ 
ables using the Cabibbo-Marinari scheme [l5j]. A sin¬ 
gle sweep is made of 1 heatbath and 3 over-relaxation 
updates of all link variables. All lattices have an in¬ 
verse temporal length 1/Lo > T c , where T c is the critical 
temperature of the theory. We have checked explicitly 


the autocorrelation times of the primary observables by 
profiting from the long Monte Carlo histories, which are 
typically made of 0(1O 5 ) sweeps. No long autocorrela¬ 
tions were observec@. For the statistical analysis we have 
blocked together the primary observables generated in 
several hundreds of consecutive sweeps, a value which is 
always much larger than the autocorrelation times mea¬ 
sured. It is important to notice that the determinations 
of Z T and z T require (see below) the computation of ex¬ 
pectation values of single local operators only. Indeed 
increasing the spatial size of the lattice does not increase 
the computational effort at fixed statistical accuracy. 

A. Determination of Z T 

The direct determination of A//A£fc in Eq. (EH1) would 
involve the computation of the ratio of two partition func¬ 
tions with different shifts at the same value of L^/a and 


2 Notice that in Ref. [l 3 | a different condition was imposed. Since 

we were interested in Z T ( 7 /q ) in a limited range of g 2 , we defined 3 At these values of Lq fluctuations of the topological charge away 
Z, r (V / 2 ) as in Eq. (12911 but at finite Lq. from zero are heavily suppressed. 
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FIG. 2. Left: the renormalization factor Z T (gg) as a function of the bare coupling for Lg/a = 3 (blue), 4 (green) and 5 
(red). Right: the renormalization factor Z T (gg) defined in the limit a/Lg —> 0 together with the fit to the formula in Eq. (1381) 
and the one-loop analytic result in Eq. H2HD- 


gg. Since the relevant phase spaces in the path inte¬ 
gral of the two systems overlap very poorly, the ratio 
cannot be estimated in a single Monte Carlo simula¬ 
tion. A possible way out is to define a series of physical 
systems with actions which interpolate between the two 
original ones, and then use the Monte Carlo procedure 
of Refs. [1;, 111, E3] ■ The calculation, however, becomes 
quickly demanding for large lattices since the numerical 
cost increases quadratically with the spatial volume. 

To bypass this problem we can profit from the fact that 
A// is a smooth function of g(j at fixed values of Lg/a 

and L/a in the range of chosen values. Its derivative with 
respect to <?q can be written as 

~dg%~K£k = 2^1 {^e-a/Lo* “ < 5 >*+«/£ofi} > ( 33 ) 

where ( S)$ stands for the expectation value of the ac¬ 
tion in Eq. (|22jl . Although the quantities on the r.h.s. 
of Eq. © have values which are close to each other, 
their difference can be computed at a few permille ac¬ 
curacy with a moderate numerical effort. The differ¬ 
ence {{S)^_ a/L(j - k - (S) €+a/Lo s} has been computed for 
£ = (1,0,0) and L/a = 48 at 63, 59 and 48 values of gg 
for Lg/a = 3,4 and 5 respectively. A sample of values is 
reported in Table |1J while all of them are shown in the 
left plot of Fig. |I] At each value of Lg/a the points are 
interpolated with a cubic spline, and the resulting curve 
is integrated over g q. The free-case value is computed 
analytically by using Eq. (IE 121) , and is added to the inte¬ 
gral. The systematics induced by the interpolation and 
the numerical integration of the data is negligible with 
respect to the statistical error. 

To complete the calculation of Z T (gg), the expectation 
value (Tyh)^ is measured in a dedicated set of simula¬ 
tions. It is computed for £ = (1,0,0) and L/a = 48 at 
66 , 60 and 38 values of gg for Lg/a = 3,4 and 5 respec¬ 
tively. A sample of values is reported in Table El and 
all of them are shown in the right plot of Fig. |I] By 
interpolating the results with a cubic spline, the renor¬ 


malization constant Z T (gg) is finally determined by the 
tree-level improved version of Eq. (EHD given by 

-MSo) = -{-XT , rib ~ free case i ■ (34) 

The resultt@ for Z T (gg) at L/a = 48 and L 0 /a = 3,4 and 
5 are shown in the left plot of Fig. [2] At the larger value 
of Lg/a = 5, discretization effects in a/Lg are within 
our statistical errors. Those due to the finiteness of a/L 
have been checked by computing Z T at L/a = 16 and 
Lg/a = 3 in the full range of gg, and at Lg/a = 5 and 6 for 
gg > 0.85. The results at Lg/a = 3 for L/a = 16 and 48 
are statistically compatible, and their central values differ 
at most by 0.5% toward the larger values of gg. Since on 
the lattices with Lg/a = 5 and L/a = 48 those effects 
are expected to be suppressed at least by an additional 
factor of 1/8, we conclude that they are well within the 
statistical errors. We thus quote the values of Z T (gg) 
at Lg/a = 5 and L/a = 48 as our best results in the 
limit a/Lg —>• 0, see right plot of Fig. [2] Even if defined 
by renormalization conditions which differ from ours by 
discretization effects, our values of Z T (gg) at gg > 0.8 
agree with those in Refs. BUI which, however, in many 
cases have a much larger statistical error. 


B. Determination of z T 

The renormalization constant z T is computed by im¬ 
posing the tree-level improved version of Eq. m given 


4 Preliminary results have been reported in Ref. [l 8 j 

5 At this small volume we have computed A//Aeither with 
the method described in this section, or with the Monte Carlo 
procedure in Ref. Q|- The numerical results are in agreement 
within statistical errors. 
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FIG. 3. Left: the renormalization factor z T (go, a/Lo) as a function of ( a/Lo ) 2 for the five values of go indicated in the legend; 
it is also shown the extrapolation to (a/Lo) = 0 with the fit function in Eq. (1371) . Right: the renormalization factor z T (g o) at 
a/Lo = 0 for the 5 values of go simulated, together with their fit to the formula in Eq. (1391) and with the one-loop analytic 
result in Eq. (|28l> . 


by 


W = 


1-4 2 

ik 






— free case 


}, (35) 


with 


L£,k 

Mi+e k ) 


5 G Z . 


(36) 


The latter condition guarantees that the WI remains 
valid at finite volume as it stands Q. The expectation 
values of and of the difference (Tpp — (T^)^ 

are measured straightforwardly in the same Monte Carlo 
simulatiorH. The free case is subtracted analytically by 
using its expression in Appendix [E] In practice we chose 
£ = ( 1 / 2 , 0 , 0 ) and q = 8 so that the ratio of the spatial 
linear size over the temporal one is fixed to be L/Lg = 20. 
We simulated 5 values of gg in the range 0 < gg < 1 with 
temporal length L 0 /a = 4,6,8 and 12. The results for 
z T (g'o ) are given in Table Hill and they are shown in the 
left plot of Fig. El Discretization effects turn out to be 
quite larger than for Z T ( g p) at the smaller values of L 0 /a. 
Our best extrapolation to a/Lo = 0 is given by the overall 
fit of the data at Lg/a = 6 , 8 and 12 to the function 


have removed the points at Lg/a = 6 from our best fit, we 
have fit each set of points independently with a quadratic 
function in (a/Lo) 2 , we have amended the combined fit 
function by adding a quadratic term in g% to the coef¬ 
ficient of (a/Lo) 2 , and we have added a quadratic term 
in (a/Lo) 2 in Eq. ( 1371 ) . The results of all these fits are 
statistically compatible with those obtained in our best 
fit to the function in Eq. (EH) and the selection of data 
points chosen. We take the maximum spread of the cen¬ 
tral values from the various fits as a systematic error due 
to the extrapolation, and we add it in quadrature to the 
statistical one. The final results are shown in the right 
plot of Fig. El 


V. RESULTS AND CONCLUSIONS 


The final results for Z T (gfy are shown in the right plot 
of Fig. El They are very well represented by the expres¬ 
sion 


^G?o 2 ) 


1 - 0.4457 gl 
1 - 0.7165 g^ 


0.2543 < 7 p 


z T (gl,a/Lo) = z T (gl) + b l gl . (37) 

The quality of the fit is very good, and it leads to the val¬ 
ues of z T (go ) shown by the black points in the same plot. 
To check for the systematics associated to the extrapo¬ 
lation, we have performed a variety of different fits: we 


It is interesting to notice that the difference (Tqq^j — (t// 

requires roughly 10 times the statistics needed for (7 ( /; to meet 
the same relative statistical error. 


+ 0.4357 g^ — 0.5221 gg (38) 

in the full range 0 < gg < 1 , a function which coincide 
with the expansion in Eq. (l28l) to order gg. The deviation 
of the curve from the data is smaller than the statistical 
accuracy, see right plot of Fig. El The error to be at¬ 
tached to Z T (gg) computed as in Eq. ( 1551 ) is 0.4% up to 
go < 0.85, while it grows linearly from 0.4% to 0.7% in 
the range 0.85 < g 2 < 1. Within our statistical errors, 
the non-perturbative determination starts to deviate sig¬ 
nificantly from the one-loop result at <?p ~ 0.25. 

Our best results for z T (gg) are shown in the right plot 
of Fig. El In the full range 0 < g'g < 1, they are well 
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represented by the expression 


Appendix B: Continuum notation 


1 - 0.5090 gl 
~ 1 - 0.4789 gl 


(39) 


In the Euclidean space-time, the path integral of the 
SU(3) Yang-Mills theory is defined as 


a function which again coincide with the expansion 
in Eq. (E51) to order g g. In this case the error to 
be attached to the values in Eq. 439) grows linearly 
from 0.15% to 0.75% in the interval 0 < < 1. 

The one loop result agrees with the non-perturbative 
determination up to g$ ~ 0.4 within our statistical errors. 

The above results for and z T (So) dearly show 

that in the range of g$ where the Wilson action is fre¬ 
quently simulated, one-loop perturbation theory is not 
adequate for computing the renormalization constants of 
the traceless components of the energy-momentum tensor 
defined as in Eq. (l26l) . Shifted boundary conditions offer 
an extremely powerful tool to compute them, and there¬ 
fore to define non-perturbatively the energy-momentum 
tensor on the lattice. The strategy implemented here can 
be easily generalized to QCD, and to (beyond Standard 
Model) QCD-likc or supersymmetric theories. 
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Appendix A: SU(3) conventions 

The Lie algebra of SU(3) may be identified with the 
linear space of all hermitian traceless 3x3 matrices. In 
the basis T a , a = 1... 8, with 


Z = 


DADcDc e s , 


(Bl) 


where the measures on gauge and ghost fields are defined 
as usual. The action is defined as 



C = C G + C GF , 


(B2) 


with 


C G = —^2 Tr 


Ffj, v 


(B3) 


C GF = ^Tr 
9o 


d„A v 


2 T 
2 ^ 
% 


d^cV^c 


where go is the bare coupling constant, Ao is the gauge¬ 
fixing parameter, the trace is over the color index and 


Fpv = d^Av - d v A il - i [A^ A v \ , (B4) 

V li c = d li c-i[A li ,c], A^ = A“T a . 

The ghost fields c and c are in the adjoint representation 
of the SU(3) group, i.e. c = c a T a and analogously for c. 


1. BRST transformations 

The action (1B2I) is invariant under the BRST transfor¬ 
mations defined as [Tiil-EiTI] 

6 Ap = 0 D^c 

Sc = A 0 9 (d^Af,) (B5) 

Sc = i6 c 2 

where 9 is an infinitesimal Grassmann constant. They 
are nilpotent up to the equations of motion of the ghost 
field c. In fact if we define 

Scj) = (j)' - <j) = 0A(j >, (B6) 


Tr[T Q ]=0, T Qt = T Q , (Al) 

the elements of the algebra are linear combinations of 
them with real coefficients. The structure constants f abc 
in the commutator relation 

[T a , T b ] = if abc T c (A2) 

are real and totally anti-symmetric in the indices if the 
normalization condition 

Ti[T a T b ] = i S ab (A3) 


where is one of the fundamental fields which transforms 
as in Eqs. (IB5I) . it is easy to prove tlratQ 

A 2 A^ = 0 , A 2 c = 0 , (B7) 

while 

A 2 c = A 0 c . (B8) 

By using Eqs. (IB7I) and (IB8I) it is easy to show that the 
BRST transformations are nilpotent, up to the equation 


is imposed. 


7 To this aim it is useful to notice that <5(D M c) = 0. 










of motion of c, when acting on any product of fundamen¬ 
tal fields at arbitrary space-time points and thus on any 
functional of them. 

The gauge-invariant part of the Yang-Mills action (IB2I) 
is BRST-invariant because the BRST correspond to in¬ 
finitesimal gauge transformations with parameter 9c(x). 
The gauge-fixing part of the action turns out to be BRST- 
invariant too. It can also be written as a BRST rotation 


of a functional plus a term which, after integrating by 
parts, is proportional to the equation of motion of c and 
serves to cancel the term coming from Eq. (IB8I) . 

2. Equations of motion 

The equations of motion for the gauge field are given 

by 



F ap \ a + d^d a A a a + -4 f abc 
% % 


where O represents a generic string of fields, and the co¬ 
variant derivative for the adjoint representation is defined 
as in Eq. (IB4I) . i.e. 

TV Fp V = dp Ffj, v - i [A^, Fp V ] . (BIO) 
The equations of motion for the ghosts are 

(l^cn*)0) = -(-^0), (BID 



Appendix C: Energy-momentum tensor in the 
continuum 


The continuum theory is invariant under the group of 
space-time translations 

x 'n= x »- e »’ <t>\ x 'p)=<t>{ x v), ( C1 ) 

where <f> indicates generically one of the fields Ap,c,c. 
The associated WIs can be derived in the usual way by 
studying the variation of the functional integral under 
local transformations parameterized by £p[x) 

SAp(x) = £ p (x) dpAp(x ), 

<5c(x) = £ p {x) dpc{x) , Sc(x) = c(x) < d p £ p (x) (C2) 
The resulting integrated WIs are 

(C3) 




SO 

6A“(x) 


(B9) 


where SO is the variation of the string of fields O un¬ 
der the transformation (1C2|) . The canonical energy- 
momentum tensor of the theory can be written as 


where 


rpC rpG,C \ rpGF,C 

[IV [AV ' fliy 5 


rpG,C _ _ rp 

n 2 
y o 

rpGF,C _ ^0 rp 
„2 
9 o 

+ 4tt 

9o 


Fp, a di/A a 


— S C G 

U I1V ) 


Aa du Ap 

(dp,c){d v c) + ( d v c)(Dp,c ) 
For e v {z) = e„S^(z — x), Eq. (IC3I) gives 


(C4) 

(C5) 

(C6) 

Sp V C FG . 


e v (dpT° v (x)0) = -{S x O) , (C7) 

and when all operators of the string O are localized far 
away from x, the classical conservation identities 


(dpT^(x)0)= 0 (C8) 

are recovered. The canonical energy-momentum tensor is 
neither symmetric nor gauge invariant. To make it both 
symmetric and gauge invariant one applies the Belinfante 
procedure, and use the equation of motion. The resulting 
tensor satisfies the on-shell WIs in Eq. (IC8I) , and it gives 
the same conserved charges of the canonical tensor when 
inserted in on-shell correlation functions. The ambiguity 
left by the use of the equations of motion allows one to 
define the energy-momentum tensor as the one derived 
by exploiting the re-parameterization transformations of 
the theory coupled to an external gravitational field [lo|, 
ITlI . |22| . All definitions related by terms which vanish 
by the equation of motion are equivalent provided the 
corresponding contact terms are taken into account in 
the WIs. The symmetric energy-momentum tensor is 
defined as 


rjiB rr\G,B _i_ rj~iGF,B 

/ns iiv ' 1 


(C9) 


d A ze v {z) (dpTp V (z) O) = — (SO) 


where 













9 


T G ] B = -* Tr 
So 


/IIS 


* [icy 


_s c G 

U UU ] 


(CIO) 


rpGF,B 

^ flU 


2 Ao 

So 


Tr 


AfidvdotAct A u d^.d a A a T d^ ^ ^d a A a dp A^ T A a d 0 


So 


Tr 

&iiC D v c H - di/C D^c 

- 5\ Tr 

d a cD a c 



So 



By comparing Eqs. (IC4I) and (IC9I) . it is quite easy to show that 


= d,T B v +8AAI 


f 1 [D a F afJ ) a + ^ d^d a A a a + \f abc (d^) c c 
9 o 9 o 


L So 


(Cll) 


r 


i.e. the two four-divergences differ by terms which are 
proportional to the equations of motion. If we insert last 
equation in the WIs (IC7I) and we use the the equations 
of motion (IB9I) we arrive to 


to T^(x) without any modification. It is worth nothing 
that the gauge-invariant energy-momentum tensor gen¬ 
erates the very same charges in on-shell correlation func¬ 
tions as all the other definitions in this appendix. 


e I/ (d„T B u (x)0) = -(6 x 0)+e v d^A“(x) J ^) . (C12) 


5A a {x), 


Appendix D: Renormalization of the action density 
in dimensional regularization 


It is also useful to notice that 


rpGF,B _ a 77 i ; — Tr 

- 1 flu — ■— 1 flU T UflU 9 11 

9 o 


ld n D n 


(C13) 


where A is the BRST variation defined in appendix | 
and 


Ifiu = — Tr 
So 


A^d v c Ajjd^c T 


(C14) 


’fyd a A a ) 


c + A a d a c ] 


When the interpolating operator O is gauge-invariant, it 
is thus appropriate to define a gauge-invariant energy- 
momentum tensor 


which satisfies 


(dfiTfiu(x) O) = (d,T B u (x) O) , (C16) 


where the term proportional to the equation of motion of 
the ghosts is null because a gauge-invariant operator is 
independent of the c field. The WIs (IC12I) applies as well 

I 


In this appendix we report the essential formulas in 
dimensional regularization which are needed in the paper, 
for a recent review see Ref. [23j and reference therein. In 
dimensional regularization one replaces f d 4 x —> f d D x, 
and renormalizes the coupling constant as 


2 _ 2e 2 ry -1 

9o — 9 9 Ag , 
where D = 4 — 2e. The /3-function is 


(Dl) 


dg 


= -eg + P(a), 


9 d 


where 


/3(s) = -g 3 ^2hg 2k 


Tfiu — T^ — 2 A^fiuF a pF a p\ (C15) and 

So k J 


bn = 


1 11 

(47t) 2 3 


■N r . 


fc=o 


bi = 


1 34 

(47r ) 4 "3 




(D2) 


(D3) 


(D4) 


with the number of colours being N c = 3 in our case. 

In presence of shifted boundary conditions it holds 

|-(Tofc) € = . (D5) 

dg L 0 dg 2 g% d£ k 

which can be written as (xq 7 ^ 0 ) 


d ,rp \ 1 dg 0 g 3 

77-\J0 fc/C = 77 - 7 ;- 

dg 2 dg g 3 0 


^(T 0k (x 0 ) {F^F^m R )z,c 


(D 6 ) 
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where Z E and the renormalized density are defined in 
Eq. m . The expectation values of the renormalized 
operators in Eq. (ID 6 I) are finite and expandable in powers 
of g. By following Ref. see also Ref. [gj, the ratio 


Since this quantity cannot have poles in e 

R = -~, (Dll) 

9 


2 e d 9 o 9 3 7 Z E Z g 

u —- Z E = — eg- -— 

d 9 9o P(e,g) 


(D7) 


must then have a series in g with no poles in e. In di¬ 
mensional regularization the coefficients of the poles in 
e in the renormalization constants start at 0(g 2 ), and 
therefore 


eZ E Z g = e + R , R = r kQ 2k > 

k =1 

which implies 

^ ZeZ g _ <?(e + R ) 
e9 p{e,g) eg~0{g) ' 

If we expand in g the denominator we obtain 
g(e + R) = (i,Rg)^(£\ k _Rg 

*9-Pis) V P ) tr'o\eg) p 


(D 8 ) 


(D9) 


(DIO) 

I 


and therefore 


eZ E Z g = e - ^ . (D12) 

9 

The renormalization constant of the energy-density oper¬ 
ator is unambiguously fixed from the one of the coupling. 


Appendix E: Lattice free theory with shifted 
boundary conditions 


In this appendix we report the results for the expecta¬ 
tion values of (T|j^)£, (T^) € - (T^) € (no sum over k), 
(7)^])£, and A//A£fc in the free theory on the latticfl 
In the infinite volume limit the expectation value of the 
momentum density is given by [ij 


8 f d 3 p sin(p 0 )sin(p fc )X] O!# 0 ifc cos 2 (p Q / 2 ) 






/ BZ 


(2tt)3 


4 sin 2 (^)+uj 2 


where 


fc =1 


t 2 |W • 2 /Pk\ 2 n£ 

= P‘£, Up = 4 2^ sin (y) , Po = -y - v P ■ 


If we notice that Q 


1 1 pix 0 (2Tv£/L 0 — <t>) 

S(^),w, x 0 ) = — Y^J 


To 4sin 2 (f F -f)+o ; 2 


1 


2 sinh Cj 


o UJX 0 


o — UX 0 


0 iLo<f>-\-LoUj _ ^ ^iLq^—LqU} _ ^ 


where w = 2 sinh {Cj/ 2), and that real and imaginary parts of E read 

sinh(Lou)/ 2 ) cosh[w(Lo /2 — £o)] — sin 2 (Lo<?i/ 2 ) sinh(Aa;o) 


Re E(<^>, w, Xq ) = 
Im E(^>, u>, Xq) = 


sinh(d)) (cosh(L 0 u>) — cos(T 0 ^)) 
— sin(Lo</>) sinh(Aa;o) 


we arrive to 


2 sinh(d>) (cosh(L 0 o)) — cos(T 0 </>)) 
d 3 p 


( T 0 k)i= 8 [ jlrT3 sin (Pk) l m S(</) p , ui p , 1) Y cos2 (y) ■ 

JBZ l 77 ) a^0,fe 


(El) 

(E2) 

(E3) 

(E4) 

(E5) 

(E 6 ) 


8 The lattice spacing is set to a = 1 in this appendix. 
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Analogously, for the traceless diagonal component of the energy-momentum tensor we obtain 

Lo — 1 j3„ 

Lo A Jbz ( 2 tt ) 3 4 sin 2 (f)+oA 


(rS) £ - (4f k = 1 d3p 


I cos(po) - cos(p fc )J A sin 2 (p a ) - cos( 2 p 0 ) - cos( 2 p fc )j A cos 2 (y)j , 

a^O ,k a^O ,k 

which by summing over l gives 

( r do>€ “ ( T kk)i = 4 / 77^{Re£(</> P ,w p , 1)^ sin 2 (p Q ) - Re £(0 P , w p , 2)^ cos 2 (y) 

^ ' c*#0,fc a^O.fc 

+ ReE(<(> p ,w p , 0 ) cos(2p fc )^ cos 2 (y) - cos(p fc ) A sin 2 (p a )J | . 

For the trace part we obtain 


a^O,k 


a^O ,k 


(E7) 


(E 8 ) 


(^ = -^E / 

0 £=0 


d 3 p Ec./^aCOS 2 ( 2 r)s in2 (P/3) 


(2tt) 3 


4sin 2 (^-) + 


which by summing over l gives 


( T d 2 J )€= 8 / |Re£(</> p ,u; p ,2)^Eos 2 (y) - Re £(0 p , w p , 1) ^ sin 2 (p fc ) 

Jbz K fc=i ^ fc=i 

-Re£(<)> p ,u; p ,0) A [sin 2 (p fc ) + cos 2 (y) + 2cos 2 (y) ^sin 2 (p g ) j . 

/C = l 

In the free theory the discrete derivative of the free energy in Eq. is given by 


At 

A£fc 

which by summing over l gives 

At 


L ° 1 ' j3 ~ w 2 +4sin 2 ((p 0 ~p k /L 0 )/ 2 ) 




t =0 




d 3 p 


In 


w p + 4sin 2 ((p 0 +Pk/L 0 )/2) 


= 2 


d 3 p 


In 


cosh(L 0 u p ) - cos (L 0 4> p + p k ) 


sz (27 r ) 3 L cosh (iow p ) - cos {L 0 4>p - Pk) _ 

i the 

E> 


All previous equations remain valid in finite volume if one makes the substitution 

f d 3 p 1 


' BZ 


(27T) 3 


V ■ 


and defines a prescription for the zero mode. 


Appendix F: Numerical results 


(E9) 


(E10) 


(Ell) 


(E12) 


(E13) 


For a representative sample of values of g q that we have in Table 

simulated we collect the results for the difference of the 
average plaquettes entering Eq. (1551) in Table HI and the 


values of ( 7(^)5 at £ = ( 1 , 0 , 0 ) in Table ITT1 The values of 
( T dfc>€ and ( T oo>€ - ( T kk)t for 4 = ( 1 / 2 , 0 , 0 ) are given 


[1] S. Caracciolo, G. Curci, P. Menotti, and A. Pelissetto, 
The energy momentum tensor for lattice gauge theories , 


Annals Phys. 197 (1990) 119. 
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p 

lg ^(ReTr U^ u }^ +a / Lo k (RcTr ^/x^)^_ a /L 0 fc] 


L 0 /a = 3 

L 0 /a = 4 

II 

53 

o 

6.0 

5.489(10) xlO -4 

3.028(7) x 10 -4 

4.2564(33) xlO” 4 

6.03 

4.886(10) xlO -4 

2.443(6) x 10 -4 

3.8484(38) xlO" 4 

6.125 

3.601(14) xlO" 4 

1.491(7) x 10~ 4 

1.0016(28) xlO" 4 

6.5 

1.576(11) xlO -4 

5.160(37) xlO" 5 

2.339(20) xlO" 5 

7.0 

7.65(8) xlO" 5 

2.232(33) xlO -5 

8.88(11) x 10 -6 

8.0 

2.96(5) xlO" 5 

7.08(25) xlO' 6 

2.43(9) xlO -6 

9.0 

1.604(38) xlO" 5 

3.60(17) xlO” 6 

1.09(11) xlO -6 

10.0 

1.041(23) xlO -5 

2.07(8) xlO" 6 

6.6(9) xlO -7 

11.0 

7.77(22) xlO -6 

1.49(8) x 10~ 6 

3.8(6) xlO -7 

12.0 

5.85(25) xlO” 6 

1.06(6) xlO -6 

3.2(5) xlO -7 

13.5 

4.29(35) xlO” 6 

7.9(8) xlO -7 

2.0(6) xlO -7 

17.0 

2.39(14) xlO -6 

3.8(5) x 10 -7 

1.3(4) xlO -7 

20.0 

1.52(12) xlO -6 

2.8(4) xlO -7 

8.2(28) xlO -8 

24.0 

1.14(8) xlO -6 

2.16(27) xlO' 7 

4.6(24) xlO -8 

30.0 

7.5(6) xlO" 7 

1.48(26) x 10 -7 

5.5(14) xlO -8 

50.0 

2.99(23) xlO” 7 

6.6(8) xlO" 8 

2.3(7) xlO" 8 

80.0 

1.25(15) xlO -7 

2.09(39) x 10 -8 

0.97(37) x 10 -8 


TABLE I. Values of the difference of the average plaquettes 
measured at bare coupling /3 = 6/go on lattices of size 48 2 3 x 
L 0 /a. 


P 



Lo/a = 3 

L 0 /a = 4 

Lo/a = 5 


(xlO -3 ) 

(xlO -3 ) 

(xlO -4 ) 

6.0 

-5.2735(27) 

-1.3772(13) 

-2.826(9) 

6.03 

-5.3921(29) 

-1.4447(11) 

-4.047(6) 

6.125 

-5.6976(29) 

-1.6064(13) 

-5.568(5) 

6.3 

-6.1359(37) 

-1.7977(12) 

-6.797(6) 

6.5 

-6.5124(28) 

-1.9495(12) 

-7.610(7) 

7.0 

-7.1554(29) 

-2.1899(20) 

-8.786(7) 

8.0 

-7.9077(37) 

-2.4488(20) 

-9.916(8) 

9.0 

-8.3673(21) 

-2.5947(30) 

-10.550(7) 

10.0 

-8.6896(14) 

-2.7002(32) 

-10.981(8) 

11.0 

-8.9385(16) 

-2.7780(31) 

-11.288(6) 

12.0 

-9.1331(20) 

-2.8358(20) 

-11.538(6) 

13.5 

-9.3654(23) 

-2.9111(16) 

-11.822(6) 

17.0 

-9.7261(22) 

-3.0181(16) 

-12.276(6) 

20.0 

-9.9253(21) 

-3.0862(17) 

-12.525(6) 

24.0 

-10.1097(15) 

-3.1420(17) 

-12.768(7) 

30.0 

-10.2941(17) 

-3.1987(18) 

-12.995(6) 

40.0 

-10.4792(18) 

-3.2587(6) 

-13.240(6) 

60.0 

-10.6608(17) 

-3.3148(7) 

-13.446(7) 


TABLE II. Values of measured at bare coupling /? = 

6/go on lattices of size 48 3 x Lo/a and £ = (1, 0, 0). 
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<?$>« - (Tj$)s 
(x 10 -3 ) 


L/a = 80 
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CD 

II 

e 
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